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Abstract: We investigate thermodynamic properties of one-dimensional U(A^) supersym- 
metric gauge theories with 4 and 8 supercharges in the planar large- limit by Monte Carlo 
calculations. Unlike the 16 supercharge case, the threshold bound state with zero energy 
is widely believed not to exist in these models. This led A.V.Smilga to conjecture that 
the internal energy decreases exponentially at low temperature instead of decreasing with 
a power law. In the 16 supercharge case, the latter behavior was predicted from the dual 
black 0-brane geometry and confirmed recently by Monte Carlo calculations. Our results for 
the models with 4 and 8 supercharges indeed support the exponential behavior, revealing 
a qualitative difference from the 16 supercharge case. 
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1. Introduction 

Large-A^ supcrsymmctric gauge theories in low dimensions are useful laboratories to test 
the AdS/CFT correspondence and its generalizations [1, 2]. Among these, the Id case 
is particularly simple since the gauge theory is just a quantum mechanical system, and 
therefore one may hope to test the duality relations explicitly and to understand them more 
in depth. Indeed Monte Carlo studies of such a system with maximal supersymmetry (16 
supercharges) have been performed by Tising non-lattice simulations [3], which reproduced 
various predictions from the dual string theory including a'-corrections in the planar large- 
N limit [4, 5, 6, 7].^ Another reason for interest in the one-dimensional case is that the 
same model, which represents the worldvolume theory of N DO-branes, has been proposed 
as a fully non-perturbative formulation of uncompactified M-theory in the light-cone frame 
[11]. Interesting observations relevant in this direction are also obtained from Monte Carlo 
calculations of correlation functions [7]. 

In this paper we apply the same Monte Carlo method to the study of non-maximally 
supersymmetric quantum mechanics with 4 and 8 supercharges.^ As opposed to the model 
with 16 supercharges, there is no known gravity dual for these models. On the other hand, 
a property common to all these supersymmetric models is that there exist flat directions 

lattice simulation of the same model has been performed in refs. [8, 9] with qualitatively consistent 
results for thermodynamical quantities. See also refs. [10] for earlier calculations based on the Gaussian 
expansion method. 

^In fact the 4 supercharge model has been studied [3, 8] prior to the 16 supercharge model for the 
purpose of testing the method. In order to address the issues given below, however, we need to study the 
system at much lower temperature, as we do in this paper. 
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in the potential, which are not hfted quantum mechanically due to supersymmetry unlike 
their bosonic counterparts. As a result, the theory contains not only the discrete states 
that the bosonic models have, but also the scattering states forming the continuous branch 
of the spectrum [12]. In the 16 supercharge model, on top of the states just mentioned,^ it 
is known that there exists a threshold bound state, which is somewhat extended in the fiat 
directions and yet has a finite norm [15, 16, 17]. Such a state, which is crucial for the M- 
theory interpretation [11], is considered not to exist in the non- maximally super symmetric 
models [15, 16, 17]. 

Recently A.V.Smilga [18] conjectured that the above difference between the 16 super- 
charge model and the other supersymmetric models may lead to a qualitative difference in 
the temperature dependence of the internal energy. In the 16 supercharge case the power- 
law behavior E oc T^^/s ^t low T was predicted from the dual black 0-brane geometry [2], 
and was confirmed by the Monte Carlo simulation [4, 6, 9]. Smilga first gave a theoretical 
argument on the gauge theory side for the particular power "14/5". Here an important 
role is played by normalizable excitations around the threshold bound state with a new 
energy scale (proportional to N~^/^) suggested from the effective Hamiltonian for the rel- 
evant O(A^) degrees of freedom in the flat directions. For the 4 and 8 supercharge cases, 
he conjectured that the internal energy decreases exponentially E oc exp(— c/T) at low T 
assuming the absence of normalizable states with a new energy scale. However, he also 
mentioned a possibility that there exist normalizable states associated with the effective 
Hamiltonian with a new energy scale proportional to N~^. In that case one obtains E ocT^ 
at low T following the same argument as in the 16 supercharge case. Our Monte Carlo 
data support the exponential behavior rather than the power-law behavior. 

As other basic properties of the supersymmetric matrix quantum mechanics, we also 
study the phase structure along the temperature axis, which turns out to be qualitatively 
the same for all the supersymmetric models. There is only one phase, in which the center 
symmetry is broken spontaneously. 

The rest of this paper is organized as follows. In section 2 we define the supersymmetric 
matrix quantum mechanics and discuss their basic properties. In section 3 we briefly review 
Smilga's argument. In section 4 we present our Monte Carlo results for the internal energy 
and compare them with the low temperature behaviors suggested by Smilga. In section 
5 we compare our results with the energy spectrum obtained in the 4 supercharge model 
for the SU(2) gauge group [19]. Section 6 is devoted to a summary and discussions. In 
appendix A we derive the expression we use to calculate the internal energy in actual Monte 
Carlo simulation. 

2. The models and their basic properties 

The supersymmetric matrix quantum mechanics are defined by the action 

S = ^ jyttr !^^{DtXif + ^i^aDt^Pa - l[Xi,Xjf - ^Va(7i)a/3[^i,V'/3]} , (2-1) 

detailed analysis of the continuum spectrum and its implications on the physics of supermembranes 
have been given in refs. [13] and [14] . 
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where Dt = dt ~ i[A{t), ■ ] represents the covariant derivative. The bosonic matrices 
A{t), Xi(t) (i = 1, 2, • • • ,d) and the fermionic matrices tpait) (a = 1, 2, • • • ,p) are N x N 
Hermitian matrices, where p = 4,8, 16 for d = 3, 5, 9, respectively. The models can be 
obtained formally by dimensionally reducing J\f = I super Yang- Mills theory in D = d+1 
dimensions to one dimension, and they can be viewed as a Id U{N) gauge theory, where 
A{t), Xi{t) and ipa{t) arc the gauge field, adjoint scalars and spinors, respectively.^ The 
pxp symmetric matrices 7^ obey the Euclidean Clifford algebra {7^, 7j} = 2Sij. The models 
are super symmetric, and the number of supercharges is given by p. 

In order to study the thermodynamics, we impose periodic and anti-periodic boundary 
conditions on the bosonic and fermionic matrices, respectively, which breaks supersym- 
metry. The extent /3 in the Euclidean time direction t then corresponds to the inverse 
temperature /3 = T~^. 

The action is invariant under the shifts 

A{t) ^ A{t) + a{t)l , Xi{t) ^ Xi{t) + Xil , (2.2) 

where a{t) is an arbitrary periodic function, and Xj is an arbitrary constant. In order to 
remove the corresponding decoupled modes, we impose the conditions 

trA{t) = 0, / dttvXi{t) = (2.3) 
Jo 

for all i = 1, 2, • • • ,d. 

Since the coupling constant g can be absorbed by rescaling the matrices and t appro- 
priately, we set the 't Hooft coupling constant X = g^N to unity without loss of generality. 
This implies that we replace the prefactor in the action (2.1) by N in what follows. In or- 
der to put the theory on a computer [3], we first fix the gauge as A{t) = ^diag(ai, • • • , ajv), 
where — vr < cti < tt, include the corresponding Fadeev-Popov determinant, and then in- 
troduce a Fourier mode cutoff A as 

A A' 
Xi{t) = ^-e--* , Mt) = E ' (2.4) 

n=— A r=—A' 

where oj = The indices n and r run over integers and half integers, respectively, and 
we set A' = A — 5 . The breaking of supersymmetry due to finite A is shown to disappear 
quickly with increasing A [3]. The fermion determinant is positive semi-definite for the 
D = A model ^ even at finite A. However, it is generally complex for the D = 6 and D = W 
models. As is done in previous works for the D = 10 model, ^ we simply omit the phase of 

*Hcrc wc use the notation of Majorana spinors to describe the D = 4,6, 10 cases in a unified manner. 
When we write a code for the -D = 4, 6 models, we rewrite the action in terms of Weyl fermions, which has 
only I components. 

^The proof goes similarly to the case of the totally reduced model [20]. 

^Previous results for the D = 10 model obtained in this way agreed well with the predictions from 
the gauge-gravity duality [4, 5, 6, 7]. It is therefore expected that the fluctuation of the phase is totally 
decorrelated with all physical quantities. While there is some evidence that this is indeed the case, complete 
understanding is still missing. 
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the fermion determinant when we study the D = Q model. Simulations has been performed 
with the Rational Hybrid Monte Carlo (RHMC) algorithm [21], which is quite standard in 
recent lattice QCD simulations. 

As a fundamental quantity in thermodynamics, the free energy T = —^\nZ{fi) is 
defined in terms of the partition function 



where the suffix of the measure [ • ]p represents the period of the field to be path-integrated. 

However, the evaluation of the partition function Z{j3) is not straightforward in Monte 
Carlo simulation. We therefore study the internal energy defined by 



which has equivalent information as the free energy, given the boundary condition = E 
at r = 0. In Appendix A we explain how we calculate the internal energy E in actual 
simulation. 

A common property of the supersymmetric matrix quantum mechanics for all D = 
4, 6, 10, which distinguishes them from the bosonic counterparts, is that they have a con- 
tinuum branch of the energy spectrum associated with the flat directions = in 
the potential. Its existence can be seen in Monte Carlo simulations as a run-away behavior, 
which can be probed by the observable 



In the high-T limit the fermions decouple, and the system becomes essentially a bosonic 
model, which is well-defined for any N > 2. The expectation value {R"^) at large but 
finite T can be obtained without any problem, and it can be nicely reproduced by the 
high temperature expansion including the subleading term [22] . As we lower T for a fixed 
N, the expectation value (i?^) decreases in accord with the high temperature expansion. 
However, at some T, we observe some cases in which the value of starts to increase 
endlessly during the simulation, which represents the aforementioned run-away behavior. 
If we use larger values of iV, we can go to lower T without seeing such a behavior. Theo- 
retical understanding of this property is provided in ref. [4] based on the one-loop effective 
action. We therefore consider that a well-defined ensemble can be defined at any finite T 
in the large- and large-A limits, and that the ensemble average corresponds to taking the 
expectation value within the Hilbert space restricted to the normalizablc states. 

In fig. 1 we plot the expectation value {R'^) obtained by Monte Carlo simulation. At 
low T and for N not extremely large, we find some cases in which (i?^) increases rapidly 
as the cutoff A is increased. This can be understood as finite A effects, which tend to lift 
the flat directions slightly [4], and hence to suppress the run-away behavior. One can also 
see from fig. 1 that larger N tends to suppress the run-away behavior. 

Our data suggest that the value of {R"^) in the above limit behaves differently for 
the D = 6 and D = 4 models. For L> = 6, it decreases monotonically (or stays almost 




(2.5) 



(2.6) 




(2.7) 



-4- 




2 2 

0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

T T 

Figure 1: The extent of space is plotted against the temperature for the D = 6 model (left) and 
for the D = A model (right). The dashed line in the right panel represents a possible log-divergent 
behavior (i?^) = a + 61og(l/T), where a = 2.71 and b — 0.566 obtained by fitting the data for 
7V = 48, A = 6, T = 0.6,-- - ,0.9. 



constant) as T decreases similarly to what is observed for D = 10 [4]. For D = 4, on the 
other hand, it starts to increase as T is lowered below T 1. This behavior is reminiscent 
of the divergence of the second moment in the totally reduced model of 4d = 1 super 
Yang- Mills theory [23, 20]. Since fermions obey anti-periodic boundary conditions in our 
Id model, the temperature T plays a role of the SUSY breaking mass parameter for the 
fermions. This gives a regularization to the second moment {B?}, which is logarithmically 
divergent with respect to the regularization parameter. Hence we obtain (i?^) ~ log^. 
Our results in fig. 1 (right) are consistent with this behavior, but more data with larger N 
around T ~ 0.4 are needed to confirm it unambiguously. 

As another property of the supersymmetric matrix quantum mechanics, let us discuss 
the phase structure along the temperature axis. For that we define the Polyakov line 

P = ^trVexp (^i J^dtA{t)^ , (2.8) 

where the symbol "P exp" represents the path-ordered exponential. It serves as the order 
parameter of the spontaneous breaking of the center symmetry. In fig. 2 we plot the 
expectation value (|-P|). In both D = 6 and D = 4 models, the results can be fitted by 
(|P|) ~ aexp^— which is a characteristic low-T behavior of a "non-confining" theory. 
This implies that the center symmetry is broken at any finite temperature similarly to the 
case of maximal super symmetry [4]. 

The bosonic models are known to undergo a "deconfining transition" at a critical 
temperature [24, 25, 26], below which the center symmetry remains unbroken.^ Therefore, 
the Eguchi-Kawai equivalence [27] (or the "volume independence" [28]) holds above the 
critical /3 in the bosonic models. As a result, the internal energy does not depend on T 
below some critical value. The physical meaning of this behavior is discussed in the next 
section. 

'^In fact there arc two phase transitions along the temperature axis as shown by Monte Carlo simulation 
[25] and by 1/D expansion [26]. 
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Figure 2: The expectation value (|P|) is plotted against the temperature for the D = 6 model 
(left) and for the D — A model (right). The solid line represents the results of the high temperature 
expansion [22] (up to the next leading order) at iV = 8 for D = 6 and at iV = 12 for = 4, and the 
dashed line represents a fit to the behavior (|P|) = aexp(— |;) characteristic to a "non-confining" 
theory. The parameters obtained by fitting are a = 1.13(1) and b = 0.215(6) for the D = 6 model 
and a = 1.03(1) and b = 0.19(1) for the D = 4 model. 



3. Brief review of Smilga's argument 

In this section we review Smilga's argument [18] for the low-T behavior of the internal 
energy. Note that the characteristic energy scale of the matrix quantum mechanics is 
Eschar ~ As, which is 0(1) in our convention A = 1 . 

Let us first start with the bosonic models, which can be obtained by omitting the 
fermions. The large- A?" behavior of such models was studied by Monte Carlo simulation in 
ref. [25]. The normalized internal energy -^E as a function of T is found to stay constant 
■^E = eq below the critical temperature. This can be viewed as a consequence of the 
Eguchi-Kawai equivalence mentioned in the previous section. Here .Evac = N'^^o can be 
identified as the vacuum energy of the bosonic model. ^ Beyond the critical temperature, 
the normalized internal energy -^E starts to grow with T. This behavior of E{T) can be 
understood from the behavior of the entropy S{E) using the relations 

Z = JdEe^^^'h-^^ (3.1) 

and E = —-^ log Z. The growth of the entropy S{E) changes qualitatively^ at some critical 
energy Ecr, where Ecr — E^ac = 0(1). Since the region E < E^ corresponds to the confined 
phase, the entropy is S{E) ~ 0(1). It is then clear that -^E = Sq below some critical 
temperature. As the energy approaches the critical value Ecr, the entropy S{E) grows 
rapidly and becomes O(iV^). (This transition has been studied in refs. [29, 30].) Since the 
region E > Ecr corresponds to the deconfined phase, it is more appropriate to describe 

*For instance, eq = 6.695(5) is obtained for the D = 10 bosonic model by Monte Carlo simulation [25]. 
Ref. [18] gives an explanation for Svac ~ 0{N^) based on the variational method. 

®Thc phase transition is found to be of second order for L* = 10 [26]. This implies that ^^^^ is continuous 
but fSi ' is discontinuous at E = Ecr- 
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the system at large N in terms of the normaUzcd entropy cr = -^S* as a function of the 
normahzed energy e = -^E. In particular, the saddle-point approximation becomes exact 
in evaluating (3.1) at A?^ = oo, and one obtains 

which gives the normalized internal energy e as a function of the temperature. This explains 
why £ has nontrivial dependence on P above some critical temperature. 

Next we discuss the case of supersymmetric matrix quantum mechanics. Here we 
consider only the normalizable states, which are included in our Monte Carlo simulation 
as we discussed below cq. (2.7). Since the supersymmetric models are not confining, it is 
more appropriate to describe the system at large in terms of the normalized entropy 
(T = -^5 as a function of the normalized energy e = ^E. 

For the D = 10 model one obtains 

£-7.4irx (3.3) 

at low T from the gauge-gravity duality, and this power-law behavior (including the coeffi- 
cient) is confirmed accurately by Monte Carlo simulation^'^ [6]. Prom eq. (3.2), one readily 
finds that 

a^e'P (3.4) 



for small e with the power 



p=-. (3.5) 
Rewriting (3.4) in terms of unnormalized variables, we obtain 

for small -^E. If we assume, for simplicity, that this formula holds even &\, E ^ 0(1), we 
obtain S ~ 0{N'^^^~'p^), which suggests that there are many normalizable states with the 
energy E ~ 0(1). (We will shortly make this argument more precise.) In ref. [18] it was 
speculated that these states should be related to the existence of the threshold bound state 
in the D = 10 model. 

Indeed one can reproduce the power (3.5) by considering the low energy excitations 
around the threshold bound state. Let us first recall that the supersymmetric models can 
have states extended in the flat directions. Wavefunctions for such states can be written 
symbolically as 

* = X{Xs\ow) '^{Xisst) , (3.7) 

using the Born-Oppenheimer approximation. Here .Xsiow parametrizes the slow oscilla- 
tions along the flat directions corresponding to the Cartan subalgebra of the gauge group, 
whereas Xfast parametrizes the fast oscillations along the orthogonal directions. Typically 



^"The power of the subleading term at low T was also derived in ref. [6] from the gravity side by considering 
the a' corrections. Monte Carlo results reproduced that power, too. 
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these states correspond to non-normalizablc states, and the spectrum becomes continuous. 
In the D = 10 case, however, there also exists a threshold bound state, which is a nor- 
malizable state with zero energy that can be expressed in the form (3.7). The excitations 
around such a state is expected to have a new energy scale E^^^ ~ N~9 [18] from the 
form of the effective Hamiltonian [31, 32] for the states (3.7). Considering that the degrees 
of freedom in the effective Hamiltonian is O(A^), the entropy should be S{E) = 0{N) for 
the energy E -^^chax ~ 0(A^9). Imposing this requirement on the power law behavior 
(3.6), one can determine the power p and finds that it agrees precisely with (3.5). 

Strictly speaking, the prediction (3.3) from the gauge-gravity duality is valid for 

_ 10 4 

N 21 ^ T ^ 1. On the other hand, here we are considering an energy scale E ~ 0(iV9), 
which corresponds to the temperature T ~ 0(iV~9) slightly below the lower bound of the 
validity region. The precise agreement therefore implies that the power law (3.3) actu- 
ally holds for a range of T wider than naively expected. Similar observation was made 
in the comparison of correlation functions with predictions from the gauge-gravity duality 
[7]. The power-law behavior of the correlation functions obtained from the gauge-gravity 
duality was found to hold in the far infrared regime, in which the supergravity analysis is 
not valid naively. 

For the D = 4, 6 cases, calculation of the Witten index suggests that the threshold 
bound state does not exist [15, 16, 17]. Therefore one may naively consider that there is 
no normalizable state written in the form (3.7). The lowest energy level^^ in the discrete 
spectrum is considered to have the energy Ey^c = 0(1), and hence we have e — > as T — > 
in the large- A?^ limit. Then the entropy is S{E) ~ 0(1) for E ~ 0(1), and it is expected to 
increase almost linearly with the energy E as it approaches S{E) ~ 0{N'^) at E ~ 0(Ar2). 
Therefore a more plausible behavior for the normalized entropy a{£) at small e is 



where the logarithmic factor is introduced to ensure e — > as T — )• 0. Then from eq. (3.2) 
we obtain 



On the other hand, the form of the effective Hamiltonian for the continuum states (3.7) 
suggests a new scale -E^^^ ~ A^^^ [33]. If there exist normalizable states with this new 
scale, the same argument^^ as the D = 10 case suggests (3.6) with p = \, and therefore 



We will see that our Monte Carlo data support (3.9) rather than (3.10). 



^^In rcf. [18] the power was determined by matching the energy at temperature T ~ E^^^ with E ~ 
NE^^^j., which is equivalent to the argument given here. 

^^In the D = 4 model with N = 2, it was shown explicitly that these states form a supermultiplet [19] 
and has a non-zero energy. See section 5 for more details. 

'^^Here it is assumed that the form (3.6) holds even at the energy scale E ~ NE^{^^^ = 0(1). 




(3.8) 



c 



(3.9) 



e ~ e T 



(3.10) 
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Figure 3: The normalized internal energy -^E is plotted against temperature for the D = 6 model 
(left) and the D = 4 model (right). The solid line represents the results of the high temperature 
expansion [22] (up to the next leading order) at iV = 8 for D = 6 and at iV = 12 for D — A. The 
triangles represent the results obtained at low T by the large-A extrapolation in fig. 4. The dashed 
lines represent the behavior (4.1) obtained by fitting the data at low T. 
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Figure 4: The normalized energy -^E is plotted against for the D = 6 model (left) and for 
the D = A model (right). The data points for each T can be fitted nicely to a straight line, which 
enables a sensible large A extrapolation. 



4. Results for the internal energy 

In fig. 3 we plot the internal energy at high temperature. The results agree nicely with the 
results of the high temperature expansion obtained in rcf. [22]. 

Next we consider the low-T behavior. Since the physical cutoff scale is given by 
uA. = 2ttTA, which is proportional to AT, one needs to increase A as one goes to lower 
T. In fig. 4 we plot the internal energy against for various T. We see that the data 
points for fixed T can be fitted nicely by a straight line. Based on this plot, we make an 
extrapolation to A = oo. 

Figure 5 shows the log plot of the normalized internal energy extrapolated to A = oo 
as a function of the temperature for the D = 6,A models. A straight line in this figure 
represents the behavior 

l^E = ae^p(^-^^ . (4.1) 
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Figure 5: The log plot of the internal energy -^E against the temperature for the _D = 6 model 
(left) and for the _D = 4 model (right). The solid line represents a fit to the exponential behavior 
(4.1) using the 3 data points with lowest T. The dashed line represents the results for the £) = 10 
model [6] obtained by fitting the Monte Carlo data to the form -^E = 7.41T"5" — cT~, where 
c = 5.58. 





Figure 6: The log-log plot of the internal energy -^E against the temperature for the D = Q 
model (left) and for the D = 4 model (right). The solid line represents a fit to the power law 
behavior (4.2) using the 3 data points with lowest T. The dashed line represents the results for the 
D = 10 model [6] obtained by fitting the Monte Carlo data to the form -^E = TAIT^ - cT^ , 
where c = 5.58. 



Fitting the data points with 3 lowest T, we obtain a = 5.5 it 1.8, b = 1.4 ± 0.1 for the 
D = 6 model, and a = 14.1 ± 0.3, b = 2.47 ± 0.01 for the L> = 4 model. The dashed hne 

represents the results for the D = W model, which cannot be fitted by a straight line at 
low T. Note also that the internal energy decreases much more rapidly in the -D = 6,4 
models than in the D = 10 model. 

In fig. 6 we show the log-log plot of the same data. A straight line in this figure 
represents 



—E = cT'^ 



(4.2) 



and the slope is given by —q. Fitting the data points with 3 lowest T, we obtain c = 3.8ibl.4, 
q = 3.5± 0.5 for the D = 6 model and c = 1.63 ± 0.09, g = 3.9 ± 0.1 for the L> = 4 model. 
The obtained values of q are much larger than 2, and we do not see any tendency that the 
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slope becomes closer to (—2) at lower T. Clearly our data are more consistent with (3.9) 
than with (3.10). 



5. Compcirison with the energy spectrum for SU(2) 



Figure 7: The histogram represents the normal- 
ized entropy cr(e) function of the normalized 
energy e in the = 4 SU(2) supersymmetric 
model [19]. The dashed line represents the result 
(5.1) obtained by fitting the Monte Carlo data 
for E{T) to the behavior (4.1). 



In ref. [19] the energy spectrum of the 
D = 4 model for the gauge group SU(2) 
was obtained by diagonalizing the Hamilto- 
nian directly. The Hilbert space was trun- 
cated by the number of bosons. The dis- 
crete states and the continuum states can be 
clearly identified by the convergence behav- 
ior as the cutoff on the number of bosons is 
increased. The energy levels for the discrete 
states were seen to converge, whereas those 
for the continuum states were seen to de- 
crease slowly without convergence. The full 
list of the discrete energy levels obtained by 
this method is given by tables 9 and 10 of 
ref. [19]. In particular, the lowest level in 
the discrete spectrum is eight-fold degener- 



On the other hand, from (4.1), we obtain 



a = Aelog {^] +Be , 



B 



log a + 1 



(5.1) 



where A = 0.405 and B = 1.48 plugging in the values of a and b obtained by fitting the 
Monte Carlo data for E{T) to the behavior (4.1). 

In order to compare the results for N = 2 with the result (5.1) for = oo, we should 
first note that the coupling constant g is set to unity in ref. [19]. The internal energy for 
arbitrary g can be obtained by E{g) = g'^/^E{g = 1). On the other hand, when we take the 
large-AT limit, we fix g'^N = 1. Therefore, we first rescale the energy spectrum obtained 
in ref. [19] by a factor of (i)^/^. Then we count the number of states (taking into account 
the degeneracy) within the interval E ^ E + SE. Dividing the number by the bin size 
5E, we obtain the density of state p{E), and hence the entropy S{E) = log p{E). The 
normalized quantities'^ a = -A^S(E) and e = -^E are defined by setting N = 2. Let us 



^^Four states form a supermultiplet, and there are two of them due to the symmetry tif 6 — rip, where 
np represents the number of fermions. 

^®We have also tried the normahzation factor (instead oi ^) motivated by the fact that the results 

in ref. [19] arc obtained for the gauge group SU(2) instead of U(2). The agreement with the large- iV result 
does not change drastically, however. 
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note that a cannot be defined as a continuous function of e at finite N. In fig. 7 we plot 
a against £ as a histogram with the bin size 5e = 0.1. At larger N, it is expected that 
the histogram becomes smoother even for smaller bin size, and that one obtains cr as a 
continuous function of e in the large- A?^ limit. The dashed line represents (5.1) for A/" = oo 
derived from our Monte Carlo results. It would be interesting to extend the method of 
ref. [19] to the N = 3 case, and to see whether the histogram comes closer to the N = oo 
curve (5.1). 

6. Summciry and discussions 

In this paper we have performed Monte Carlo studies of supersymmetric matrix quantum 
mechanics with 4 and 8 supercharges at finite temperature. Similarly to the 16 supercharge 
case studied previously, the potential for the bosonic matrices has a flat direction, which is 
not lifted quantum mechanically due to supersymmetry. This is seen in our Monte Carlo 
simulation as the run-away behavior. However, it is possible to suppress this behavior 
and to define a well-defined ensemble by using sufficiently large N. We consider that 
the ensemble obtained by our simulation corresponds to restricting the Hilbert space to 
normalizable states. 

While the phase diagram turns out to be similar to the 16 supercharge case, we observe 
a notable difference in the behavior of the internal energy. It decreases much faster as the 
temperature decreases than in the 16 supercharge case. In fact, our Monte Carlo data for 
the internal energy are consistent with exponential decrease at low temperature in striking 
contrast to the power-law behavior in the 16 supercharge case. This behavior was predicted 
by Smilga assuming the absence of normalizable states with a new energy scale unlike in 
the D = 10 model. Thus our results provide independent evidence for the peculiarity of 
the D = 10 model suggested previously in the literature [15, 16, 17]. 

To our knowledge, there are no concrete proposals for a gravity dual of the non- 
maximally supersymmetric models studied in the present work. We hope that our explicit 
results in the planar large- A" limit would provide a useful guide for constructing such an 
example. That would give us new insights into the gauge-gravity duality in non-maximally 
supersymmetric cases. 
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A. Derivation of the formula for the internal energy 



In this Appendix we explain how we calculate the internal energy defined by (2.6). In 
ref. [22] it is shown that the internal energy E is related to the expectation values 

^ = (fb) + (^f) , (A.l) 
where the operators and £f are defined as 

^^ = -\^ j^dttv (V'alTOa/sI^i, V'/j]) • (A.3) 

However, the calculation of {£{) is time-consuming since it requires the construction of 
the fermion matrix and inverting it. In fact we can use an identity to trade it off with a 
quantity written solely in terms of bosonic variables. This trick has been first proposed in 
ref. [8] and it was used also in refs. [4, 6]. 
Let us consider the change of variables 

Xi^e'Xi{t) (A.4) 

in the path integral representation of the partition function (2.5). One finds that the 
path-integral measure transforms as^^ 

and each term in the action transforms according to the order of Xi. Since this is just a 
change of variables, the value of the partition function should not change. This implies 
that the 0(e) terms should cancel. Thus we obtain an identity 

^ i^jyt^^ \^^\{DtXif - ^\[Xi,Xjf - ^Va(7i)a/3[^i, V'/j]}^ = d{iV2(2A + 1) - 1} . 

^ ° (A.6) 
Note that the derivative term on the left-hand side has linear divergence at A — >■ oo, which 
is given by the 0(A) term on the right-hand side. Using this identity, we obtain 



F- ^ 



(5b)-^d{(2A + l)iV2-l} 



(A.7) 



The advantage of using this formula instead of (A.l) is that the evaluation of Sb does not 
take much time. The disadvantage is that the fluctuation of S'b is 0(\/A), and therefore 

one needs to increase the statistics linearly in A to keep the statistical error constant. This 
is related to the fact that the relation (A.7) involves the cancellation of the 0(A) terms. 
We do not encounter any problem, however, in the parameter regime investigated in this 
work. 



^The "—1" in the exponent is due to the constraint (2.3) 
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